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Abstract. Sticltjes integral theorem is more commonly known by the phrase 'integration by 
parts' and enables rearrangement of an otherwise intractable integral to a more amenable form; 
often permitting completion of an integral in closed form. Applying Stieltjes integral theorem to 
the Hilbert transform yields an alternate integral definition, which is homeomorphic and exhibits 
accelerated computation. By virtue of the convolution theorem, the integral is mapped to Fourier 
image space and delineates requirements for the inverse Fourier transform, also, these requirements 
reveal the underlying reason for increased computational speed. Lastly, an alternative to Cauchy's 
integral formula is deduced by extending the line integral with logarithmic kernel into the complex 
domain. 



Key words. Stieltjes integral theorem, Hilbert transform, integrals with logarithmic kernels, 
Cauchy's integral formula, Kramers-Kronig, causal analytic signals, Harmonic theory. 

o . 

^vq ; AMS subject classifications. 42A38, 42B10, 42E05 



1. Propaedeutics. The Hilbert transform is not uniquely defined in the litera- 
ture and differs by a minus sign; to avoid confusion, the Hilbert transform is defined 
in Definition (jl.ip and was adapted from R. N. Bracewell, also, the descriptor 'of the 
first form' will become evident later in this essay [1]. 

Definition 1.1 (Hilbert transform of the first form). The classical definition 
for the Hilbert transform (Ji?) of a function f(x) on domain T> is defined to be the 
following definite integral: 

where, the Cauchy principal value V is taken at all discontinuities. 

The Hilbert kernel, (x — a;') -1 , experiences a discontinuity for every case where 
| x' = x; this property classifies the Hilbert transform as an improper integral |S]. The 

function f(x) may possess additional undefined points within the domain; regardless, 
being an improper integral bars treating the Hilbert transform as a Riemann inte- 
gral. For an arbitrary function f{x), integration often necessitates expanding into the 
complex domain; whereupon, Cauchy's integral theorem proves indispensable, thus, 
k^J ■ V represents the Cauchy principal value. 

For an n-dimensional Euclidean space, the Hilbert transform can be applied along 
each dimension, T> C 1", provided the integrals exist. The Hilbert transform can also 
be applied to complex domains, T> C C™, once again, assuming integral convergence. 
If a finite number of poles exist with no accumulation points in the complex domain 
for complex function f(z) E C", then the function is said to be meromorphic [2]. A 
closed Jordan curve removing all poles from the domain forms a piecewise holomorphic 
function that represents the original meromorphic function 2 . Holomorphic functions 
are continuous, analytic functions. Cauchy's integral theorem may be applied to 
holomorphic functions to evaluate the Hilbert transform of the function f(z) in the 
upper half of the complex plane. 
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Consider a real- valued constant function f(x) over the interval a < x < b, the 
resulting integral is null, see equation (|1.2p . If interval (0,S) were deducted from 
the interval of integration for the Hilbert kernel, the value of the integral is -co, no 
matter how small S is made j8] . Direct numerical integration of an arbitrary function 
f(x) must judiciously handle the symmetric requirements in the neighborhood of 
any singularity within the domain; because of this behavior, the Hilbert transform is 
unstable and difficult to integrate numerically [7]. 



The Hilbert transform as a convolution integral is shown in equation (|1.3|) . where 
an asterisk (*) represents convolution. The Fourier transform of the Hilbert 
kernel is the sign function, written as sgn(s); in addition, the function is multiplied 
by the imaginary number, i = \J — 1. Fourier theory states the convolution of two 
functions is equivalent to multiplying the Fourier image of each function in Fourier 
image space. Taking advantage of the convolution theorem, the Hilbert transform is 
mapped to Fourier image space as the multiplication of the sign function times the 
Fourier image of the function f(x), finally, the inverse Fourier transform (j?"^ 1 ) is 
applied to return to the original domain [I] . 



The Fourier image of the Hilbert transform replaces the numerical instability of 
the Hilbert kernel by the unbounded sign function. Depending upon the nature of 
the function f{x), Shannon's sampling theorem may or may not be easily satisfied 
[B] . Because the sign function is unbounded, all discrete transforms risk numeric error 
associated with under-sampling. 

The function f(x) and the corresponding Hilbert transform of f(x) comprise a 
pair of functions referred to in unison as a Hilbert pair. In combination, the Hilbert 
pair form what is generally known as an 'analytic signal' or 'strong analytic signal'; 
the pair of functions are heavily relied upon in signal processing and causal physical 
theory, e.g. electromagnetic theory supports the well-known Kramers-Kronig relations 
[3]. Each function of the Hilbert pair represents an orthogonal geometric trajectory 
and in tandem form a set of functions completely covering a Riemann manifold [5] . 

2. Assembly. It is common practice to apply, then immediately reverse an oper- 
ation to find a new mathematical identity; such is Stieltjes integral theorem. Stieltjes 
integral theorem is more commonly known as 'integration by parts' and is stated in 
Theorem ([2.1)1 . which has been adapted from D. V. Widder's codification [8]. Integra- 
tion by parts is frequently employed to transform an otherwise intractable integral to 
a form that is more amenable to known techniques of integration. For the purposes of 
this essay, let the first and second instances of integration in Stieltjes integral theorem 
be referred to as "the first form" and "the second form", respectively. 

Theorem 2.1 (Stieltjes integral theorem). If integrand function f(x) and inte- 
grator function a(x) are continuous on the interval (a, b), then the sum of integration 
for the first and second form equals the limit sum of the product f{x)a{x). 




(1.2) 




(1.3) 




(2.1) 
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Proof: Take the derivative of the product f(x)a(x) immediately followed by 
integration, i.e. J d(f (x)a(x)) . Applying the product rule for derivatives yields two 
terms, whereupon integration distributes and yields the sum of two integrals, finally, 
the limit sum is a consequence of the Fundamental Theorem of Calculus [8]. 

The Hilbert transform, Definition (jl.ljl , can be identified as the first integral form 
in Stieltjes integral formula. To find the Hilbert transform of the second form, we first 
identify the integrator function a(x) to be the indefinite integral of the Hilbert kernel, 
which is the natural logarithmic function and leads to the following intermediary form: 

^{f{x)} = ^ y f(x')d\n(x - x') 

Reversing the role the integrand and integrator functions play yields the Hilbert 
transform of the second form, see Definition (|2.2|) , where the convolution form of the 
integral is also shown. The definition consist of the derivative of the function f{x) 
multiplied by a shifted logarithmic kernel, then integrated over the entire real line. 

Definition 2.2 (Hilbert transform of the second form). The Hilbert transform 
of the second form (Jt?(2)) is defined by the following definite integral over the entire 
real line: 

J? {2) {f(x)} = ^ J°° \n{x-x')^ldx' = Un{x)*df{x) (2.2) 

where, the Cauchy principal value V is taken at all discontinuities. 

An interesting property of the convolution theorem arises when conjoined with 
the derivative theorem, namely, if one attempts to take the derivative of the entire 
convolution integral, one is free to operate the derivative on one function singly and 
then convolve with the other function, which has been untouched by any previous op- 
eration; the choice of which function to operate the derivative is completely arbitrary. 

Theorem 2.3 (Derivative of a Convolution Integral). The derivative of a con- 
volution is equal to the convolution of either function with the derivative of the other 

(/*<?)' = /'*<? = /*</ (2.3) 

where prime signifies a first order derivative. 

Theorem (|2.3[) provides an alternate proof for the Hilbert transform of the second 
form. Beginning with the convolution form of the Hilbert transform of the first form, 
the derivative is extracted, then applied to the function f(x), also, the derivative is 
signified by prime rather than the italicized letter d. The omission of a minus sign 
is intentional; retention of the minus sign leads to the negative of the result shown 
below. The absence of the minus sign indicates the Hilbert transform of the second 
form is, in fact, equal to the negative of the Hilbert transform of the first form. 

— * fix) = - ( Hx) * f(xj) ' = - Hx) * f{x) (2.4) 

7TX TT \ J TT 

The Hilbert transform of the second form can be mapped to the product of two 
Fourier image functions in Fourier space by invoking the convolution theorem. It is 
required to know the Fourier transform of the natural logarithmic function to com- 
plete this transformation and this affords another opportunity to implement Stieltjes 
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integral theorem. The symmetric Fourier transform of the natural logarithmic func- 
tion is stated on the lefthand side of equation (|2.5I) . Employing Theorem (|2.ip . the 
integral involving the logarithm and Fourier kernel can be rearranged and equated to 
the righthand side of equation (|2.5p . 



1 r°° 1 1 

ln(x)e- 27rixs dx = / -e~ 2 ™ s dx ln(x)e- 27rl:ES 

2nis /_~, x 2nis 



(2.5) 



The integral on the righthand side of equation (|2.5I) transforms 1/x to — i7rsgn(s). 
The limit sum of the product between the natural logarithmic function and the Fourier 
kernel would appear, at first glance, ambiguous; because, the Fourier kernel oscillates 
and does not converge at infinity. The ambiguity can be managed by considering 
either the limit inferior or limit superior of the functions involved. 

Evaluating the limit sum gives rise to two terms: the first term is evaluated 
approaching the infinite point along the positive direction of the real line; the second 
approaches minus infinity from the negative direction. In the case of evaluating the 
logarithmic function at minus infinity, the principal value of the natural logarithm 
is introduced, yielding the absolute and argument (arg) of the independent variable; 
finally, substitution leads to the second and third terms found on the righthand side 
of equation (|2.6[) . The limit superior (lim) of both the first and second terms cancel. 
The limit superior of the last term, the argument (arg) times the Fourier kernel, is 
equal to iir. 



ln(or)e 



= lim hn(x)e- 27rzxs ) - lim ( In |or|e- 27r<xs )- (2.6) 

x— ►oo I J x— y— oo I J 

lim { i axg{x)e~ 27Tlxs \ => -i?r 

X— > — oo I J 



The sign of the limit involving the argument multiplied by the Fourier kernel is 
ambiguous, for the limit inferior would return —iir; hence, the solution to the limit 
sum is, in general, ±in. Collecting results, the Fourier transform of the natural 
logarithmic function can be displayed succinctly in equation (|2.7j) . where the double 
arrow signifies a simplification of the algebraic formula after realizing the following 
equivalent relationship: sgn(s)/s = l/|s|. 
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Using equation (|2.7|) . the Fourier image of the Hilbert transform of the second 
form involves the Fourier transform of the natural logarithmic function and the Fourier 
image of the first derivative of the function f(x), stated in equation (|2.8[) . Since the 
Fourier image of the derivative of a function is equal to 2iri times the Fourier image of 
that function f(x), we can swiftly dispense with the troubling prospect of taking the 
derivative of what could be a discrete set meant to represent the continuous function 
f{x). Multiplication by 2itis and subsequent algebraic simplification leads to the last 
equality specifically, right of the double arrow in the following: 



{*&>{/(*)}} = ^ (r 4) ± ^ 



As the transform variable s varies over all real numbers, {s\s £ M}, the algebraic 
term s/\s\ + 1 takes on the value of 2 for all positive real numbers or zero for all 
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negative real numbers; this can be compactly represented by twice Heaviside's step 
function, written as 2H(s). Heaviside's step function evaluates to 1/2 at the origin 
and the algebraic term evaluates to twice unity at the origin, hence, we must add 
a suitable term to account for the additional value missing, namely, the product 
between the imaginary number, Dirac Delta function S(s) and the Fourier transform 
of the function f(x). This allows further simplification of equation (|2.8|) and leads to 
equation (|2.9|) . 

!( - 2iH(s) - i<5(s)W{/(a;)},for s/\s\ + 1 
(2.9) 
(2iH(-s) + i<5(s)W{/(x)},for s/\s\ - 1 

The Hilbert transform of the second form replaces the sign function with, essen- 
tially, twice Heaviside's step function; effectively nullifying either the entire negative 
or positive real line, depending upon which form is employed. As a consequence, an 
alternate perception of the Hilbert transform is that the transform suppresses either 
all negative or positive frequencies of a signal; thus, either advancing or retarding the 
phase of a signal. 

The inverse Fourier transform operated on both cases found in equation (I2.9P 
would retrieve the Hilbert transform of the second form in the original domain. To 
that end, consider that the inverse Fourier transform of Heaviside's step function 
equals the sum of the Hilbert kernel and the Dirac Delta function S(x): 

l v ') 2ttx 2 

Applying the above identity for Heaviside's step function to both cases found 
in equation (|2.9[) produces the results found in equation (|2.10j) . The initial result 
states that the function f(x) is to be convolved with the Hilbert kernel with either 
the addition or subtraction of that function convolved with Dirac's Delta function, 
finally, either addition or subtraction of the imaginary number times the function f{x) 
itself. 

jq 2 ){f(x)} = — * f(x) ± iS(x) * f(x) ± if(x) (2.10) 
K ' l > irx 

Further simplification is attained by realizing the convolution of any function with 
the Delta function equals that function alone, namely, 5(x) * f(x) = f(x). Therefore, 
the inverse Fourier transform of equation (|2.9p can be reduced to equation (|2.11j) . The 
earlier indication that the Hilbert transform of the second form was somehow related 
to the Hilbert transform of the first form can now be stated explicitly as a Corollary. 

Corollary 2.4 (Equivalence of the first and second form). The Hilbert trans- 
form of the second form is equal to the negative of the classical Hilbert transform 
definition with either the addition or subtraction of function f(x) multiplied by 2i. 

{/(*)} = -^{f(x)} ± 2if(x) (2.11) 

where the argument of the logarithmic function is responsible for the ambiguity in 
sign. 

Substituting Corollary (|2.4p for the definition of the Hilbert transform of the 
second form in Steiljtes integral formula, Theorem (|2.ip . reduces the relation to a 
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statement concerning the infinite point. Evidently, the residue at the infinite point 
is equal to 2mf(x); because of the ambiguity of the logarithmic function, the sign of 
the product is either positive or negative. 

COROLLARY 2.5 (Principal value) . Substitution of Corollary \2.J$ reduces Stielt- 
jes integral formula for the Hilbert transform to the following statement: 

^{f(x)}+J? {2) {f(x)} = 

je{f(x)}-je{f(x)} ± 2if(x) = i ln(z - x')f(x') 



ln(x - x')f(x') 



= ±2mf(x) 



(2.12) 



Finally, using Corollary \2A\ and Definition (I2.2p . it is possible to deduce an 
alternate form for Cauchy 's integral formula, which states the line integral around a 
shifted logarithmic kernel and the first derivative of a function f(z) is equal to the 
function f(z) without differentiation. 

Lemma 2.6 (Line integral with logarithmic kernel). Let T be a closed Jordan 
curve in domain T>, then the line integral of the derivative of function f(z) with a 
shifted logarithmic kernel is equal to the function f{z): 



(2.13) 



An n-th differential form of equation (j2.13[) is possible, much like Cauchy 's integral 
formula; albeit, this is less interesting, for it only states the integral of f^ n+1 '(z) is 
equal to f^ n '(z), thus, erasing only one degree of differentiation. In contrast, the 
derivative formula for Cauchy 's integral formula finds far more utility. It is further 
possible to apply successive derivatives to the logarithm, but this only returns a 
formula similar to Cauchy's integral formula; thus, providing little difference overall. 

3. Exhibition. Real-world signals are inherently real-valued and more than 
likely represented by discrete sets. Because discrete sets are usually encountered, 
it is most natural to implement Fourier transforms to emulate the Hilbert transform. 
Accelerated computation of the inverse discrete Fourier transform is immediately ob- 
vious for a list of elements one half the length required by any other means or schema. 
Since the Hilbert transform of the second form nullifies one half of the real line, then 
the prospect of accelerated computation of the Hilbert transform exists. 

Inspecting Table (I3.1[) . there is an immediate computational savings from halving 
the length of any given list of data points with regard to Fourier inversion; but, any 
apparent benefit swiftly diminishes as the length of the list increases. This can be ex- 
plained by the fact that regardless of the algorithm used the difference in the number 
of operations required for the inversion becomes negligible to the total required num- 
ber of operations as elements are increased; nevertheless, measured statistics clearly 
show a marked increase in all cases tested. 

Lists possessing a length equal to an integer power of two, generally speaking, fall 
under the class of algorithms referred to as the Fast Fourier Transform. As evidenced 
by the data in Table (|3.1I) , it is apparent that lists of length not equal to an integer 
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Table 3.1 

Statistics for computational runtimes involving both forms of the Discrete Hilbert Transform 
(DHT), where the first form is the classical Hilbert transform and the second form is Definition 
12. 2\ ). The comparative percent increase for the 2nd-DHT algorithm is relative to that of the lst- 
DHT algorithm. For an indicated trial number, the average runtime to complete the inverse discrete 
Fourier transform (DFT) is shown with accompanying standard deviation. The power N determines 
the number of elements operated on by the forward DFT; the length of the list is equal to N — 1 for 
the inverse discrete Fourier transformation in the case of the 2nd-DHT algorithm. 



Integral 
form 


Power 
2 N 


Trial 
number 


Percent 
increase 


Mean runtime 
(milliseconds) 


Standard Deviation 
(milliseconds) 




10 


3000 




2.367 


5.600 


Ist-DHT 


12 


3000 




9.204 


104.800 




12.5 


300 




3.269 


0.770 




18.5 


300 




569.509 


42.930 




20 


300 




3283.840 


185.500 




10 


3000 


515 


0.385 


2.561 


2nd-DHT 


11 


3000 


582 


1.472 


4.562 




12.5 


300 


1364 


0.223 


0.417 




18.5 


300 


1616 


33.180 


6.893 




20 


300 


347 


740.877 


60.865 



power of two took less time to compute; an explanation of this behavior is impossible 
without being privy to internal proprietary knowledge of Mathematica's algorithms. 

Because the internal workings of Mathematica 's InverseFourier function is not 
privy, it was decided to use the average of several successive runs to remove propri- 
etary optimization methods from biasing any of the tests. Be it as it may, the percent 
increase in computational speed is relative and not dependent upon the specific work- 
ings of the algorithm employed nor computer language, hardware architecture, &c. 
Regardless of the specific runtime for each length of list entertained, the relative 
computational increase in runtime was faster for the second Discrete Hilbert Trans- 
form algorithm (2nd-DHT) over the classical Discrete Hilbert Transform algorithm 
(lst-DHT). 

The infinity norm of the difference between numeric output from both algorithms 
were within machine error and determined both numerically identical, specifically, 
| |e| |oo < 15.9546-logarithm base ten scale; furthermore, this was true for any test 
function employed. Runtimes were calculated using Mathematica's built-in Timing 
function, intended for this purpose. All numeric calculations were written and exe- 
cuted using Mathematica version 5.2 software, from Wolfram Research, installed on 
an AMD Phenom™ Triple-Core Processor 2.30 GHz machine, with 2.00 GB RAM 
memory. 

4. Epilog. The Hilbert transform enjoys prominence in Harmonic theory and 
sundry uses spanning many disparate disciplines in both theoretical and engineering 
applications. With the aid of Stieltjes integral theorem, a new novel integral form is 
found for the Hilbert transform, dubbed 'the Hilbert transform of the second form' 
for a lack of a better term. The second integral form involves an integral over the 
entire real line with a logarithmic kernel, of which, this form of integral often arises 
in many engineering and physical applications [9l [3] . 

Applying the Hilbert transform generates the Harmonic conjugate of causal real- 
valued signals. In unison, the Hilbert pair constitutes an analytic signal and embodies 
all potential information contained in the measurement, i.e. amplitude, phase and 
frequency. Many well-behaved continuous functions are susceptible to closed form 
transforms; albeit, such signals are ideal and serve mostly theory. In general, measured 
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signals are represented by a discrete set of data points; furthermore, it is to this sort 
of data representation that the Hilbert transform is more likely brought to bear and 
the underlying continuous functions rarely gleamed. 

It happens that many numerical methods exist to evaluate the Hilbert transform; 
yet, given the discrete nature, it is natural and most efficient to implement Fourier 
transforms in exercise. Benchmark tests show computational savings for the Discrete 
Hilbert Transform of the second form; furthermore, savings are solely due to the 
length of the list being truncated before the inverse Fourier transform is applied. 

Finally, by extending the domain to include a closed Jordan curve within the 
complex domain, the formula for the Hilbert transform of the second form may be 
cast as a contour integral; thus, providing an equivalent formula for Cauchy's integral 
formula, but in this case involving a logarithmic kernel. The utility of an alternate 
Cauchy integral formula becomes evident for integrals with logarithmic kernels; ei- 
ther enabling direct integration or providing a transform to an integral explicitly not 
involving a logarithm [9l [3] . 
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